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Abstract 



O . The note presents a classification of the relevant distinct types of solutions 

I of the general Friedmann equation without assuming a priori restrictions for 

If^ • the parameters occurring in this equation. The emphasis is on the case of a 

non- vanishing cosmological constant. The classification uses algebraic criteria. 
I The result is: There are four distinct basic types of models. Explicit formulas 

■ for decisive terms are given. Characteristic mutual relations of cosmological 

constant, mass and radiation density to distinguish between the models are 
calculated. 



1 Introduction 



The Friedmann equation is a model for the description of a universe satisfying Ein- 
steins field equations of the general relativity theory and Einsteins cocmological prin- 
(Nl I ciple, expressed by the ansatz of a Robertson- Walker metric. The present paper 

describes the distinct types of solutions of this equation, using algebraic criteria for 

■ the classification. Explicit formulas for decisive terms are presented. We do not 
assume a priori restrictions for the parameters occurring in this equation. 

For convenience a brief outline of the derivation of the Friedmann equation from 
the mentioned starting point is given (see e.g. [1]). Einsteins field equations read 

■ {Rfi,u - i:gtJi,yR) + J^Qi^^u + i^T^^u = 0, K := — - 



c3 ■ where {g^^u} denotes the metric tensor, {R^^u} the " Verjiingung" of Riemanns curva- 

ture tensor, R := g'^^'^R^^^, where the sum convention is used, {T^^u} the energy tensor 
and A the cosmological constant. G is the gravitational constant and c the velocity 
of light. Further Einsteins cosmological principle says: The universe is homogeneous 
and isotropic w.r.t. space. The ansatz 

ds^ = c^de - Ritf { + r\de'' + sin^ Qdcj)'')] (2) 
y 1 — er^ J 

of Robert son- Walker for the metric tensor corresponds to this principle, r, 6', (j) are 
dimension-less and time-independent coordinates. The dynamics is determined by the 
scale-factor R{t) only. It has the dimension of a length, e is the curvature parameter 
with values e = 0,+!,— 1 corresponding to the euclidean, spherical and hyperbolic 
case, respectively. The energy tensor is given by 

Tf,,y := (p + ^) u^Uu - g^,,uP, (3) 
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where p, P, Un denote mass density, pressure and 4- velocity, respectively. The cosmo- 
logical principle requires that p and P are homogeneous w.r.t. space, i.e. 

p^p{t), P^P{t), ii = {c, 0,0,0} 

where t is the eigen-time. According to (2) and (3) one obtains 

{T^,.} = diag |pc^ PPV^ PB?PsinH^ . (4) 

Putting (4) and (2) into (1) one obtains for the 00-component of the field equations 

-^-AP=--^(pc^ + 3P)P (5) 

and for the space-components 

RcPR 2 (cLrV ^ , „ 47rG/ 2 r.i 

+ ^Hr +2e-Ap2=^(pc2-P)p2. (6) 



c2 cit2 c'^ydt ) 



c 



From (5) and (6) one obtains 



Idp 3 dR / P\ 

• (7) 



c Pc V 

2 

Using the approximations P = and P = ^ for the " non-relativistic approximation 
of the present universe" and a "radiation-dominated universe", respectively, then 
from (7) one obtains 

2 

Pmat{t)R{tf ^: Amat for P = 0, p,t,(i)P(i)^ =: forP=^, (8) 

where Pmat, Pstr dcnotc mass and radiation density, respectively and where Amat and 
Astr are constants. Assuming that there is no coupling between radiation and matter 
one puts 

P Pmat + Pstr- (9) 

Then from (5) and (6) one obtains the Priedmann equation 
where 

SttG , SttG , , , 

a:^—^A^at, cr —^Astr, a > 0, a > 0, e = 0, ±1. (11) 

Further restrictions of the parameters A, e, a, a are not used at the moment. 

A is considered as a fundamental natural constant connected with the basic struc- 
ture of space-time. Since the time of the appearance of the Einstein-de Sitter model 
1932 (see [2]), where A = is assumed, the case A = was generally preferred. 
However recently the case A > was again taken into account because of new ob- 
servations. Therefore the focus of this note is directed to this case. First a brief 
overview on the case A = is given. 
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2 The case A = 



In this case the Friedmann equation (10) reads 

-ec^ e = 0,±l (12) 




2.1 The case e = 0,-1, a + > 

We put e = -6, 6 = 0, 1. Then 

aR + a + Sc^R^ > 0, R>0, 

i.e. there is no restriction for R w.r.t. possible solutions. In particular, there is no 
constant (stationary) solution R{t) = const. Either always R'{t) > or R'{t) < 0. 
In the first case R{-) is strongly monotonically increasing and it turns out that it 
increases asymptotically for t ^ oo like const ■ t for e — —1 and like const ■ t^l^ for 
e = 0. The second case R!{t) < can be excluded. 

The case e = is called the Einstein- de Sitter -model. If one introduces the 
so-called Hubble-parameter 

then, using (11), (8), (9), the equation (12) can be written in the form 

Since the present value Hito) is approximatively known, one can define a so-called 
critical density 

_ 2,H{toY 

2.2 The case e = +1, a + o- > 

In this case (10) reads 

= ■^{aR + a-c'^R'). 
R"^ 

If a = then there is an exceptional radius Rq := The region R > Rq is & 

forbidden region for solutions because in this case a — (?R^ < 0. The straight line 
R{t) — Rq represents an unstable stationary solution. The region < i? < i?o is 
allowed. There are two different types of solutions, the first one with R'{t) > starts 
at i? = and finished at i? = Rq, the second one reversal. That is, each point of 
the constant solution R{t) = i?o is a branching point for solutions, refiecting the 
unstability of this stationary solution. 

If q; > 0, then again there is an exceptional radius. 





with the same implications as before. 
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2.3 Criterion for the value of the curvature parameter e 

Writing the equation (12) in the form 

H{tf = -^P{t) - (13) 
the following criterion is obvious: 

e equals +1,-1,0 iff p{to) > Pcr , p(^o) < Per , p{to) = Per , respectively. (14) 

3 The case A > 

Obviously, there is a large distinction between the cases A = and A > because 
A is the coefficient of the "highest term" in the right hand side of (10). Therefore it 
would be difficult to start with a perturbation theoretic ansatz. For convenience we 
put 

P := -Ac", 7 := ec^ 



such that (10) now reads 



We put 



That is 



where 



p{R) - r^-1e' + ^r+^. 



a. 

p{R) = Rq{R) + -, 



As already mentioned, A is considered as an so far unknown fundamental constant 
and e is also a fixed (but also a priori unknown) constant. Thus the solutions of 
(15) depend on a and a. Obviously a classification of them depends strongly on the 
properties of the polynomial p of the fourth degree. Note that the "shape" of p is 
independent of a. This fact suggests to use the polynomial q and the parameter a as 
preliminary parameters for the classification. 

The strategy for the classification is to start with the discriminant D of q, given 

by 

(ay 27a'P - 47^ 

The parameter D is decisive for the properties of the roots of q. Therefore, first we 
recall the following facts from the theory of zeros of polynomials of the third degree, 
applied to q. 

(i) li D < then q has two different real simple roots < i?i < i?2- Note that 
a > implies Ri > 0, but for a = one obtains Ri = 0. 
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(ii) li D — then q has a positive double-root R1 — R2. 

(iii) If D > then q has no non-negative real root. 

Note that (16) implies that in the cases e = and e = — 1 always D > follows. This 
means that for these cases always the case (iii) is true. This implies that also the 
polynomial p has no positive real roots for all cr > 0. This suggests first to separate 
this special case. 



3.1 The case e = 0, -1 and a > 

In this case D > follows and consequently there are no restrictions for > 0, i.e. 
there are no forbidden regions or separating straight lines. 



3.2 The case e = +1 and = 

In the case e — one has to take into account all three possibihties for D. This 
suggests to start with the case a = 0. 

The roots R2 in the case (i) can be calculated explicitly. The result is 

where 

/ \ -3/2 . 

where 

TT 

— l<cos0<O or TT > > — . 

The limit case = 7r/2 corresponds to the case a = and the limit case = tt to 
the case = 0, i.e. to the case of maximal a. In the latter case, where Ri = i?2, one 
has cos f = I such that in this case 

R,^R2^ A-V2. 
Note that in this case, where € — +1, one has 



^mat 



M 
2^' 



(19) 



where M is the total mass (dust) of the universe, because in this case the volume of 
the 3-sphere is 27r^i?(t)^. Therefore, 

cos(/)= -^MA-^/' (20) 



and D < means 
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This means that in the case o" = 0, D < 0, one has the following situation: The 
interval R G -R2) is a forbidden region for solutions because in this case p{R) < 0, 
i.e. there is an upper region R > R2 and a lower region < R < Ri, where solutions 
are possible, except for the case a — 0: in this case the lower region disappears. 
Moreover, in the case a > one obtains two different unstable stationary solutions 
R{t) := Ri and R{t) := R2. 

A lower region appeared already in the case e = +l,A = 0, q;>0 (see 2.2). Recall 
that for a = the exceptional radius is Rq = ac~^. The new feature in the case 
A > consists in the appearance of a completely new upper region where solutions 
arc possible. If wc consider the limit A ^ then R2 ^ 00 and Ri ac~^. This 
follows from (17) and (18). One has 

(f - f ) _ iRi 

— cos 3a 
If ^ ^ then (18) yields ^ 7r/2. Now 

<j>^TT/2 — COS (j) 3 

that is I = 3^1im^^o-Ri- 

The special case a = 0, together with A > 0, cr = as before, is a limit case in 
so far as in this case the "new" upper region is present but the "old" lower region 
disappeared. In this limit case (10) reads - with e = +1 

The straight line 

separates the upper region R > Rq from the forbidden region R < Rq. Again one 
has the unstable stationary solution R{t) := Rq. The solutions in the upper region 
represent expanding universes where for i — > 00 the solution R(t) is asymptotically 
given by 

exp 

These solutions are due to W. de Sitter (see [3, 4]). They are called de Sitter universes. 

In the case a = 0, D = the forbidden region degenerates to a straight line, i.e. 
the point Ri — R2 is a 3-fold branching point for solutions because an initial point 
{^0, Ri = R2} has three possibilities to evolve: to remain stationary, to enter the 
upper or the lower region. 

If (T = and D > then there is no restriction for > 0. This means that the 
solutions can be divided into two disjoint classes, defined by R'{t) > or R'{t) < 
for all t. The latter class can be excluded. 
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3.3 The case e = +1 and a >0 

First let D < 0. Then one obtains immediately from the structure of p that if a 
increases then the forbidden region decreases, i.e. 

(Ti < 0-2 implies Ri{ai) < Ri{a2) < R2{(^2) < R2{(^i) 
This means: There is a critical value Ucr such that 

-^1(0"^) = R^ic^cr) ='■ Rcr- 

The critical radius Rcr can be calculated explicitly. This term depends only on the 
derivation of the polynomial 

Again this polynomial has two different positive real roots < it!^ < R2 which realize 
the maximum and the minimum of the polynomial Rq(R), respectively. Further 
Rcr — -^2- The calculation of Rcr yields 



Rcr ^2 1-^] cos ^ = (y) ' cos %, (21) 



1/2 

where cos^" = :^cos0, i.e. ^ > ^ f- This means, in the massless case one has 
= |, as before, and in the limit case of maximal mass = or'0/3 = 7r/4, which 
implies that in the limit case one has 



R\ — R2 — R, 



cr 



Then a simple calculation yields 

Rcrq{Rcr) = COS - ■ COS — , 

and the critical cr- value is acr '■— —^RcrQiRcr), i-e. 

(7cr = — :— cos — -COS — . (22) 
A 3 3 ^ ^ 

That is, if < then one has to distinguish between three cases: 

(i) a < acr- Then the interval (i?i((T), i?2(c)) is a forbidden region, 

(ii) a — acr' The forbidden region degenerates to a straight line, 

(iii) a > acr- There is no restriction for R > 0. 
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4 Summary 



The foregoing discussion of several cases can be summarized as follows. For simplicity 
the branching cases D — 0, a — and D < 0, a — are omitted, also the case 
a = and the cases where there is no time interval with R'{t) > 0. There are four 
distinct types of solutions: 

(A) e = or e = -1. Further a >0. (Recall D >0 in this case.) 

(B) e = +1 and: L'>0, C7>0 or D<0, ct> acr- 

(C) e — +1 and: D < 0, cr < acr- (Solution in the upper region.) 

(D) € — +1 and: D < 0, a < acr- (Solution in the lower region.) 

Next as the counterpart the essential analytic-geometric characteristic properties of 
the cases (A)-(D) are described. We choose to = as the starting time. 

(A) R{0) = 0, for t — 7> oo the radius (scalar-factor) R{t) behaves approximatively 
like exp(v^t), there is a turning point for R{-); H{t) > is strongly mono- 
tonically decreasing. 

(B) -R(O) = 0, for t — )■ oo the radius R(t) behaves approximatively like exp(v^t), 
there is a turning point for i?(-); H{t) > 0, there is a minimum point tmin 
for H{-), H{t) > H{tmin) for ^ > tmin and there is a turning point t^^H > tmin 
for hI). 

(C) R{0) = R2{o'), for t — )■ oo the radius R{t) behaves approximatively like 
exp(v^t), -R'(O) = 0, H{t) > is strongly monotonically increasing. 

(D) R(0) = 0, there is a maximum point tmax for R{'), Ritmax) = -Ri(c); H{t) > 
is strongly monotonically decreasing. 



5 Conclusions 

A celebrated observation on the large-scale structure of the universe was the discovery 
that the universe expands, the scale-factor R increases. This is an implication of 
the observation that the Hubble-constant (the Hubble-parameter of the presence) is 
positive. 

Recent observations on very far distant Supernovae, e.g. the Supemovae Cosmol- 
ogy Project (SCP) or the observations by WMPA-sateUites, show that the expansion 
is even accelerating (see e.g. S. Perlmutter [5]). This means that H{-) is neither 
constant nor decreasing at present, i.e. in a former epoch H was smaller than to-day. 
This led - in contrast to the long accepted Einstein-de Sitter-model (see [2]) - to the 
conclusion that one has to assume rather that A > 0. This assumption is compatible 
with the interpretation of the expansion as an intrinsic property of "space", which 
leads consequently to an exponential growth or increase of it (see e.g. Lemaitre [6]). 

If these mentioned observations will be confirmed by further observation projects, 
then case (A) could be excluded, i.e. then necessarily e = only a spherical model 
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could be the right one. Also (D) could be excluded, i.e. the so-called BigBang- 
BigCrunch solution. 

On the contrary, the models (B) and (C) are compatible with the mentioned 
observations. 

The model (B) represents a BigBang solution, whereas the model (C) can be con- 
sidered as an Anti-BigBang solution. Interestingly enough, this model is essentially 
the model of Lemaitrc (see [6]), it does not have a singular origin but it starts as a 
static Einstein- Universe (see e.g. [7]) and can be described for t ^ oo asymptotically 
by the de Sitter-Universe (see [3,4], see also [8]). 

That is, these observations lead to the alternative between the BigBang-model 
(B) and the Anti-BigBang model (C). At present the BigBang- hypothesis is widely 
accepted, which is supported by the 3K-background radiation and the distribution of 
Helium and Hydrogen in the universe. However, still it remains a hypothesis. 

If it would turn out by further observations in the future, for example that in far 
former epochs the Hubble parameter H was much larger than today, i.e. that H was 
decreasing in that time, then this would be a proof for the BigBang-model (of course 
within the framework of the Friedmann equation), because in this case H must have 
passed through a minimum and only in the model (B) there appears such a minimum 
for H. 

According to the Summary the alternative between the models (B) and (C) can 
be expressed by their characteristic parameters as 

D > 0,OT D < and a > Ucr versus < and u < acr- 

Finally it should be mentioned that also in the case A > the Friedmann equation 
(10) can be formally written in the form (13) if one puts 

Ac^ 

P(t) := Pmat(t) + Pstr(t) + — ^. 

OTTCt 

Then again the alternative (14) remains true. However, since all models with e = 
0,-1 are now excluded one can state 

p{to) > Per = „ ^ , ^0 present time. (23) 

ottG 

According the foregoing analysis (23) does not imply that the BigBang-BigCrunch 
model (D) is true. 



6 The special case D > 0, a = 

In this special (B)-model it is easy to calculate the characteristic points = 
R(t^),Rmin = R{tmin),H{tminT,Rw,H = R{tw,H) and H{tw^HY (see [9]). First note 
that 

Rw ^ Rmin ^ Rw,H 



because oi D > 0. Using the parameters M, A, c, G one obtains 

2GM 



7^ . / TT C 

D > means A > 
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Further 



TT C 



2 \ 2GM „ 3GM 



-^tu — ' ^ 2 GM I ' ^^^^ — ^ ^2 ' ^w,H 



TT C"^ TT C"^ 



3 I \2GMJ r ' ' 3\ 21\2GM) 

7 Some calculations for the Anti-BigBang model 

In this case one has D < and a < Ocr- We are interested in the case of nearly 
maximal mass M. According to (20) one has 



MAV2 = |^(-cos(/)), |<0<7r. 



For convenience we put 



2^ ^ 

^ := 2 ~ ^' X > 0, X small. 



Then f = f — f ^^"^ = T ~ ^ hence cosV' = ;^(— cos ^ + sin Further one 
has COS0 = v^cosV', i.e. 

Since % is small, in the first approximation w.r.t. x, i.e. neglecting higher terms than 
X, one has 

MAV^ = f^(l-|). (24) 

Further one has cos(| — f ) = -^(cos | + sin |), hence cos^(| — f ) ^ ~ x) ^-^d 

2 V' 2V^ 1 1 
cos - • COS — = - smx(l + smx) = 2 ■ X 

follows in the first approximation. Thus one obtains 

(7cr = • X- (25) 

The result is: If (24) is satisfied where x > is small and a < c^A~^X: then the 
corresponding Friedmann-model is an Anti-BigBang model (C). 

Finally, we calculate the minimal radius in this case for the limit case a = a^r- 
According to (21) one gets 

"'"VaJ ^°^3"VAj cos(--x) = A /(cos- + sm-), 
i.e. in the first order approximation one has 



VA V 2 
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Using (11) one obtains from (25) 



A 



str,cr 



47r2 [2G 




) 



X 
A' 



In the limit case D = 0, i.e. x = 0) this means 



MA'/' = ^ 
2G 



and in this case one obtains 




The constant A* |^ 



is connected with k (cf. equation (1)) by /i ■ k — Air'. 
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